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ABSTRACT 

The survey phase of the Kepler Mission includes a number of hot subdwarf B (sdB) stars to 
search for nonradial pulsations. We present our analysis of two sdB stars that are found to be 
g-mode pulsators of the V1093 Her class. These two stars also display the distinct irradiation 
effect typical of sdB stars with a close M-dwarf companion with orbital periods of less than 
half a day. Because the orbital period is so short, the stars should be in synchronous rotation, 
and if so, the rotation period should imprint itself on the multiplet structure of the pulsations. 
However, we do not find clear evidence for such rotational splitting. Though the stars do show 
some frequency spacings that are consistent with synchronous rotation, they also display mul- 
tiplets with splittings that are much smaller. Longer-duration time series photometry will be 
needed to determine if those small splittings are in fact rotational splitting, or caused by slow 
amplitude or phase modulation. Further data should also improve the signal-to-noise, perhaps 
revealing lower amplitude periodicities that could confirm the expectation of synchronous ro- 
tation. The pulsation periods seen in these stars show period spacings that are suggestive of 
high-overtone g— mode pulsations. 
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1 INTRODUCTION 

Subdwarf B (sdB) stars are horizontal-branch stars with masses 
~0.47Mq, thin (< 10~^Mq) hydrogen surfac e layers, and 
effective tempera tures from 22 000 to 40 000 K jHebeJ 1 19841 : 
ISaffer et"al]|l994b . They represent the aftermath of the core he- 
lium flash, and as such are important stars for understanding this 
rapid phase of stellar evolution. Some of these stars are nonradial 
pulsators, and we can exploit this by using the tools of asteroseis- 
mology to probe their interior structure, as well as determine ac- 
curate estimates of their total mass, shell mass, luminosity, internal 



rotation, and the extent of radiative levitation and gravitational dif- 
fusion. 

Pulsating sdB stars come in two varieti es. The first type to b e 
found was the short-period V361 Hya stars l lKilkennv et al.|[l997h . 
They are typically referred to as sdBV stars. Their pulsations are 
p-modes, with periods that range from two to four minutes, but 
they can range up to ten minutes; pulsation am plitudes are typ- 
ically near 1 per cent of their mean brightness ( iKilkennvllioOTi : 
iReed et ^l2007bl) . The second vari ety are longer-period ff— mode 
pulsa tors, first discovered in 2003 jOreen et ^ l2003l : iReed et al.1 
|2004 and known as V1093 Her (or PG 1716) stars. They are com- 
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monly referred to as long-period sdBV (IpsdBV) stars. Their pe- 
riods range from 45 min to more than 2 hr with amplitudes typi- 
cally < 0. 1 per cent though amplitudes up to 0.5 per cent have been 
observed 10stensen 2009). They also are all multimode pulsators. 
Examples of the large number of pulsation modes determined us- 
ing ground-based observations are PG 1338-1-481, with 13 period- 
icities jRandall et al. 2006b), and PG 1627-1-017 with 23 periodic- 
ities I IRandall et al.i i2006a). There is also a gr oup of hybrid stars 
that pulsate with both short and long periods dOreiro et alj|2004l : 
ISchuh et al.ll2006l) . 

Those stars that are p-mode pulsators allow us to probe the 
outer layers of these stars to below the transition from hydrogen 
to helium. For understanding their cores - and directly probing 
the regions where helium is u ndergoing fusion - we need accu- 
rate measurements of g-modes dCharpinet et al.ll200(]h . The obser- 
vational challenges in accurately identifying the pulsation modes 
are significant. This usually requires extensive photometric cam- 
paigns, preferably at several sites spaced in longitude to reduce 
day /night aliasing. However, Kepler provides long-duration contin- 
uous, homogenous, evenly spaced time-series photometry, making 
it an ideal instrument for asteroseismology. 

The striking advantage of Kepler data is the many improve- 
ments over ground-based data, particularly for long-period sdBV 
variability. From the ground, one can only observe a few pulsa- 
tion cycles during nighttime hours. If multisite observations are 
obtained, then differing instrument sensitivities become an issue, 
and because of weather, gaps will doubtlessly appear in the data. 
Additionally, ground-based observations have to contend with at- 
mospheric transparency variations that can often act on the same 
time scales as the pulsations in the stars. In contrast, Kepler data 
are optimal for these types of pulsators in that nearly gap-free data 
are obtained at a roughly constant cadence and with no atmospheric 
issues. The combined advantage in signal-to-noise and coverage 
allows us to detect many more pulsation frequencies than Earth- 
bound telescopes can obtain, even with significant effort. 

Because of their asteroseismic potential, sdB stars (and white 
dwarfs) are targets that form a key component in the Kepler 
Mission's asteroseismic investigation. The Kepler Mission science 
goals, mission design, an d overall performance are reviewed by 
iBorucki et al.l (2010) and iKoch et a l. (2010). Asteroseismologi- 
cal analysis for the Kepler Mission is being conducted through 
the Kepler Asterosei smic Science Consortium (KASCfl - see 
iGilliland et all ( 1201 Oah . Compact pulsators, including sdB stars and 
white dwarfs, are the responsibility of the Compact Pulsators Work- 
ing Group (WGll), which is overseeing the analysis of the sur- 
vey phase data on pulsating hot subdwarfs. WGl 1 target selection, 
spectroscopic properties of the targ ets, and overall results of the 
search for pulsators are described in l0stensen et alj ( l2010al . Paper 
I). The first sh ort period sdBV star found by Kepler is presented in 
iKawaler et al.l (i2O10, Paper II). Long-period sdB puls ators found 
in the first part of the survey phase are described in iReed et al] 
1 2010l Paper III) with d etailed modelling of one of those stars in 
Ivan Grootel et al.l ( 1201 Ol Paper IV). 

One of the lingering mysteries in stellar astrophysics is the 
mechanism by which sdB stars form. There are many of these 
stars, all of which have very thin surface hydrogen layers. Their 
relatively high effective temperatures signify surface hydrogen lay- 



^ At the time of this work, the Kepler Asteroseismic Investigation, which 
manages the KASC, was led by R. Gilliland, T. Brown, J. Christensen- 
Dalsgaard. and H. Kjeldsen. 



ers ranging from 10~ to 10~^Mq. One of the leading forma- 
tion scenarios relies on binary m ass transfer as the sdB progen- 
itor climbs the first giant branch dHan et al]|2002l |2003|) and, in 
fact, many sdB stars (p ulsators and non-pul s ators) show evidence 
of close companions jMaxted et al 200 ll : iMorales-Rueda et"al] 



l2003l : iNapowtzki et alj|2004 iHebed (20091) . In several cases, pul 



sating sdB stars are in close binaries with orbita l periods of a 
few hours to 1.2 days or so ("Kilkennv et al. 1998; Billeres et al.l 
l200d ; lO'Toole et alj|2004l : [Oeier et al. 2010) . Because the orbital 
period is so short, the stars should be in sync hronous rotation 
( IZahnlll975L[l977l ; lGoldreich & Nicholsonlll989h . Thus these sys- 
tems can help test ideas of tidal synchron ization through spec - 
troscopy (" Geier et al .120 id) and photometry dBloemen et al.l2O10h . 
Asteroseismic studies sug gest that the sd B component of the bi- 
narie s Feige 48 (Van Gro otel et al J I2OO8I P^9h) and PG 1336- 
018 dCharpinet et al .1120081. P = 2.4h) are indeed in synchronous 
rotation. In Feige 48 , the companion is probably a white dwarf 
dO'Toole et ai]|2004l), vi^hile i n PG 1336-018 the companion is an 
M5 star dVuckovic et al]|2009t) . 

In this paper we examine two long-period sdB pulsators 
discovered during the survey phase of the Kepler Mission: 
KIC 02991403 and KIC 11179657. Both of these stars show evi- 
dence for low-amplitude variation revealing a close binary com- 
panion with an orbital period of 9-10 hr. The photometric varia- 
tions suggest a lower-mass companion displaying the reflection ef- 
fect. The parameters of the binary systems (beyond the orbital pe- 
riod and amplitude of photometric variation) are at the moment not 
known, so in this paper we concentrate on analysis of the pulsa- 
tions and defer a detailed discussion of the binary properties to a 
later paper. We note that another complex pulsating sdB star in a 
close binary was uncovere d in the survey p hase, with preliminary 
analysis presented in .Ostensen et al.l d2010bh . 



2 OBSERVATIONS 

This paper describes data obtained by Kepler during 2009. The data 
released to the KASC Compact Stars Working Group (WGll) are 
short cadence (SC) data with an average cadence of 58.85 s. This 
pipeline provides time series of "raw" flux values for each target, 
and "corrected" fluxes based on preliminary estimates for contam- 
ination by nearby stars. Because the contamination correction will 
be improved with further upgrades to the data pipeline, in this paper 
we work with the raw fluxes provided to the KASC. For analysis of 
the flux variations, we consider fractional variations away from the 
mean flux (i.e., differential intensity A///). Amplitudes are given 
as milli-modulation amplitudes (mma), with lOmma correspond- 
ing to 1.0 per cent. 

While the data from Kepler are nearly continuous, there are 
occasional gaps in the data as the result of safing events and other 
brief anomalies. The window function of the power spectrum is 
essentially a sinc'^ function with a width equal to the inverse of the 
run length, since the duty cycle of the reduced light curves is so 
high, and the outlying points are scattered more-or-less randomly 
through the light curves. The time series data provided from the 
SC data frequently show instrumental artefacts that are at multiples 
of the long-cadence (LC) readout rate of 1/30 the SC rate. This 
produces peaks in the temporal spectrum at n ■ 566.44 /xHz, and 
experience shows that the larg est-amplitude artefac t peaks are in 
the 4000 to 7000 fiKz region tGilliland et alj|2010bl) . All variation 
frequencies in these data are well short of 566.4 /xHz and so we are 
not influenced by this artefact. 
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The two targets were selected using the process described in 
Paper I; further details about these two targets can be found in that 
paper. These two stars are the faintest of the eight pulsating sdB 
stars uncovered in the first half of the survey phase. 

KIC 02991403 This target was selected based on photometry 
from the SDSS/SEGUE survey tStoughton et al. 20021: lYannv et al.l 
l2009h . With a Kepler magnitude Kp of 17.136, this star is a spec- 
troscopically confirmed sdB star with T^s = 27, 300 ± 200 K and 
log(? = 5.43 ± 0.03. The spectroscopic parameters of this star, 
and those of KIC 11179657 quoted below, were estimated from 
spectra obtained with the ISIS spectrograph on the William Her- 
schel Telescope, having a resolution of 7?~ 1600. For details on 
the spectroscopic determination of these quantities, see Paper I. 
This star lies within the long-period sdB instability region in the 
H-R diagram. Time-series data were obtained during the first quar- 
ter year of Kepler science operation (Ql) over a 34 d span between 
BJD 2454964 and 2454997 (2009 May 12 and 2009 June 14). This 
provides a formal frequency resolution (1/T) of 0.35/iHz. Since 
Kepler data prior to Q2 were delivered without a barycentric cor- 
rection, we applied the BJD correction using the documented pro- 
cedures provided by the Kepler program 0. The Kepler Input Cat- 
alog (KIC) lists a contamination factor of 0.601 for this star. As 
mentioned earlier, we used the raw flux for KIC 02991403. When 
contamination effects are applied, the amplitude will be higher than 
the raw amplitude by a factor of (1 — c)^^, where c is the contam- 
ination factor. The KIC estimate of c is expected to be improved 
upon in subsequent analyses and data releases. The raw data con- 
tain a significant number of outlier points - we removed these points 
by first boxcar filtering the data to remove the orbital effects, and 
then identifying, for later removal, points that were in excess of 4 
times the RMS deviation of points from the mean. This removed 
1177 points from the raw light curve (roughly 2 per cent of the to- 
tal number of data points), leaving 47,939 data points for analysis. 
Finally, a slow decrease in overall flux was present in the Ql data 
on this star, so we fitted and removed a 2nd-order polynomial trend 
from the data. 

KIC 11179657 We identified this star as a potential target based 
on p hotometry from the AIS survey of the GALEX satellite 
([Martin et al. 2005). This star has a Kepler magnitude Kp of 
17.065. Spectroscopic determination (Paper I) provides an effec- 
tive temperature of 26,000 ± 800 K and log 5 = 5.14 ± 0.13. 
This places it within the long-period sdBV region of the H-R dia- 
gram, but at the low-gravity end. Kepler data were obtained during 
the second quarter during a 27 d span between BJD 2455064 and 
2455091 (2009 August 20 to 2009 September 16). This provides 
a 1/T frequency resolution of 0.43 /iHz. A total of 39,810 science 
images were scheduled. However because of safing events, loss of 
fine guidance, and other small glitches 1,492 images were either not 
obtained or were contaminated in some way. As for KIC 02991403, 
we de-trended the slow variation in the mean level of the data, and 
removed 4-a outliers, resulting in a total of 38,272 data points in 
the reduced light curve. The KIC lists a contamination factor of 
0.129, thus the amplitudes we determine in the time-series analysis 
using the raw fluxes are closer to the true amplitude in the star than 
for KIC 02991403. 



^ J. Van Cleve (ed.) Kepler Data Release Notes 2, available at: 
|http://archive.stsci.edu/kepler/data-release.html| 



3 ANALYSIS 

Our analysis proceeded in a straightforward way as has been use d 
in prior papers in this series jKawaler et alj201ol : lReed et alj2010h . 
We calculated a temporal spectrum (a Fourier transform; FT) which 
provides initial estimates of pulsation frequencies and amplitudes. 
We simultaneously fit and prewhitened the data for each of these 
frequencies using non-linear least-squares (NLLS), beginning with 
the highest amplitude frequency until we reached the 4-(t detection 
limit. We determined a by finding the mean value of the temporal 
spectrum in the frequency range of interest, excluding obvious peri- 
odicities. We begin here with a discussion of the binary light curves, 
and the assumed impact that rotation should have on the pulsations 
as we search for evidence of rotational splitting at the orbital / co- 
rotation frequency. We follow with analysis of the higher-frequency 
pulsational periodicities in the two stars. 



3.1 Binary light curves 

Both of these stars show variations with periods of about 10 h. 
KIC 02991403 and, to a lesser extent, KIC 11179657, show equal 
amplitude, peaked maxima and shallow minima, typical of the ir- 
radiation effect produced when an sdB star heats up the face of an 
M-dwarf orbiting the subdwarf (with its rotation period likely syn- 
chronized to the orbital period). These effects produce light curves 
with a dominant peak at the orbital frequency, and a smaller peak at 
the first harmonic. The bottom panels of Figure[T]show the binned 
light curves for the two stars. The specifics of the binary properties, 
including spectroscopic observations, will appear in a subsequent 
publication. 

No eclipses are evident in either light curve, and the ampli- 
tude is quite low compared to other known sdB-l-dM binaries. For 
comparison, the ampli tude of the reflection effect in th e sdB-l-dM 
binary HE 0230-4323 ('Ko enll2007l : rKilkennv et alj201 0^ is about 4 
per cent in V, and the orbital period is 10.8 h. The eclipsing sdb-l-dM 
binary PG 1336-018 ( orbital period of 2.4 hr) is almost 10 per cent 
jKilkennvet ai]|l998l) . Since PG 1336-018 is an eclipsing system, 
the inclination must be close to 90° and so represents the intrin- 
sic reflection effect amplitude. On the other hand, if the amplitude 
of the effect drops with the square of the orbital separation, then 
the intrinsic amplitude of the reflection effect is approximately 6.3 
times smaller, or 16 mma. If that is the case, then these systems may 
be closer to edge-on, since the amplitude of the reflection effect that 
we see is then comparable to what we would expect. In the oppo- 
site extreme, if the intrinsic amplitude is similar for these somewhat 
longer-period binaries, then the low amplitudes we see (down by a 
factor of 10) could indicate a much smaller inclination of 5° or so. 
From the information we currently have about these stars, it is diffi- 
cult to say what, if any, constraints the binary light curve places on 
the orbital inclination. Time-series spectroscopy will be required to 
constrain the companion type and inclination angle further. 

In the left panels of Figure|2]and Figure[3]we show the tempo- 
ral spectrum of the frequency region containing the binary period 
and first harmonic. The right panels in these figures show the area 
of the temporal spectrum where pulsations appear. 

For this analysis, we are primarily interested in how the bina- 
rity impacts the variations attributed to pulsations. For close bina- 
ries such as these, one expects that the stars will be in synchronous 
rotation - that is, the rotation frequencies of the stars should match 
the orbital frequency of the binary (il = /orb, where Q. is the ro- 
tation frequency and /orb is the orbital frequency), and that the 
rotation and orbital axes are aligned. For nonradial pulsation, rota- 
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Figure 1. The light curves of two sdB binary stars. For each star, the top panel shows an 8-day sainple of the SC 
light curve. In the bottom panels, the time-series data (with the pulsations filtered out) are folded at the indicated 
binary period, and binned into 300 s phase bins. Error bars represent the scatter about the mean for each phase 
bin and include, in part, the pulsations. The left panel shows data for KIC 02991403, and the right panels show 
KIC 1 1 179657. 



tion should cause a shift in the observed oscillation frequencies for 
modes with nonzero values of the azimuthal quantum number m. 
In particular, for slow rotation (Q, <^ u, where v is the pulsation 
frequency) it is well known that: 

Vn^l^m ^ yn,l,Vl + mQ.{l — Cn^l) (1) 

( lLedoux|[T95lh . Here, the indices n, I, and m correspond to the 
radial order n and the indices of the spherical harmonic functions 
Yj™. The quantity Cn,i describes the contribution to the splitting 
caused by the Coriolis force, and is a function of the eigenfunction 
for the appropriate nonradial mode. Rotation, therefore, can break 
the degenerate m frequencies into their 2i + 1 components sepa- 
rated approximately by the rotation frequency. As such, multiplets 
can be used to associate pulsation modes (described by n, Z, and 
m) with observed frequencies. Seismic models of stars provide the 
theoretical oscillation frequency spectrum, which can be quite rich; 
two similar models can show similar frequencies, but with different 
values for n, I, and m. Therefore mode identifications are needed 
to constrain stellar models, which in turn are used to infer internal 
conditions of the stars themselves. 

For the long-period pulsations seen in these stars, we expect 
that they are high-order g-modes. A useful approximation for these 
modes, where the horizontal motions are significantly larger than 
the radial motions, is that 

Consulting evolutionary models in the rele vant mass, TeS and log .q 
ranges from lKawaler & Hostlej ( l2005h and lCharpinet et alj bOOA . 
we find that C„,i ranges from 0.465 to 0.496 for I = 1 gr-modes in 
the observed period range (see below) and between 0.158 to 0.165 



for / = 2 modes. So the above approximation appears reasonable 
for these stars. We note that if there is some tidal deformation of the 
sdB stars in these binaries, that there may be additional frequency 
shifts associated with that change in stellar shape (i.e. one would 
expect equally spaced multiplets for all modes with the same spac- 
ing if the tidal axis coincides with the pulsation axis). For this initial 
analysis, we assume that the pulsators are spherically symmetric. 

While rotation can lift the m degeneracy and produce 21 + 1 
peaks for each value of n and I, the viewing geometry (as well 
as the intrinsic pulsation amplitude) will alter the observed am- 
plitudes. The geometric effects, at leas t, can be com puted in a 
straightforward way (see, for example. |Pesnelll l l985h . For these 
stars, which show orbital effects in the light curve, we expect that 
the stars are in synchronous rotation at the orbital period. Thus, 
/ — 1 modes are likely to appear as multiplets split by 0.51 to 0.54 
/orb with up to 3 peaks for each overtone n. For / = 2 we should 
look for a multiplet splitting of about 0.84/orb, with up to 5 peaks 
for each overtone. 

3.2 KIC 02991403 

The temporal spectrum of KIC 02991403, prior to prewhitening, is 
shown in Figure[2l in the top panels. There are several artefact fre- 
quencies in the 4500-5500 /iHz range, but most of the peaks in the 
FT of this star occur in the 100-400 /iHz range. We determine a 4- 
a detection limit using apparently variation-free frequency regions 
above and below this range. Between 1 and 150 /iHz, the 4-a limit 
is 0.40 mma, and between 400 and 550 fiHz we find a value for 
the 4-a limit of 0.31 mma. This detection limit is shown as a blue 
(mostly horizontal) line in Figure|2] 

The middle and bottom panels of Figure |2] show the detec- 
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Table 1. Frequencies, periods, and average amplitudes for KIC 02991403. 
The binary period and first harmonic are included as well. Formal least- 
squares eiTors are in parentheses in this and subsequent tables. 



ID Frequency [/iHz] Period [s] Amplitude [mma] 



Binary Period and first harmonic 



/orb 


26.1240(0.0011) 


38278.9 (1.6) 


17.55 (0.10) 


2 /orb 


52.2466 (0.0084) 


19140.0(3.1) 


2.33 (0.10) 




Pulsation Frequencies 




fl 


157.101 (0.025) 


6365.35 (1.01) 


0.484 (0.063) 


f2 


194.675 (0.028) 


5136.77 (0.75) 


0.556 (0.064) 


f3 


195.144 (0.031) 


5124.41 (0.81) 


0.513 (0.064) 


f4 


230.552 (0.022) 


4337.42 (0.42) 


0.553 (0.063) 


f5 


231.115 (0.030) 


4326.84 (0.55) 


0.416 (0.065) 


f6 


264.437 (0.012) 


3781.62 (0.18) 


0.980 (0.065) 


f7 


270.739 (0.032) 


3693.59 (0.43) 


0.382 (0.066) 


f8 


284.149 (0.020) 


3519.28 (0.25) 


0.601 (0.063) 


f9 


285.309 (0.019) 


3504.97 (0.23) 


0.661 (0.064) 


flO 


296.303 (0.036) 


3374.93 (0.41) 


0.332 (0.063) 


fll 


308.188 (0.028) 


3244.77 (0.30) 


0.434 (0.063) 


fl2 


309.268 (0.016) 


3233.44(0.16) 


0.783 (0.065) 


fl3 


334.233 (0.030) 


2991.93 (0.27) 


0.450 (0.063) 


fl4 


334.818 (0.013) 


2986.70 (0.11) 


1.144 (0.063) 


fl5 


335.390 (0.026) 


2981.60(0.23) 


0.544 (0.064) 


fl6 


369.013 (0.030) 


2709.93 (0.22) 


0.402 (0.067) 



tion limit and the prewhitening process for KIC2991403. The high- 
est amplitude variation has a frequency near 26.1 /jHz and a pe- 
riod near 0.44 d. This frequency and its associated overtone are 
attributed to binarity and shown in the left panels of Figure [2] 
We detected 16 other frequencies that we attribute to long-period 
gr— mode pulsations. These are listed in Table [T] along with their 
amplitudes and periods and indicated in Figure|2]with arrows. No 
peaks remain in the FT above the detection limit after removal of 
these periodicities. 

Table [T] lists the resulting periodicities. Of the 16 pulsation 
frequencies in the table, there are 5 sets of closely spaced multi- 
plets. We list these close spacings, along with the spacings expected 
for synchronous (and uniform) rotation in Table |2] The smallest of 
these spacings, 0.47 /^Hz, is a bit larger than the formal frequency 
resolution of 0.35 /iHz. The other close spacings have values of ap- 
proximately 0.55 to 0.58 /iHz or multiples thereof. If these small 
spacings are interpreted as a rotational splitting of I = 1 modes, 
then the implied rotation period (using Equation 1 and 2) is 10.4 d. 
If they are I — 2 sectoral m = ±2 modes, then the rotation period 
is 34 d. In either case, if this splitting is indeed caused by rotation, 
then this star is far from synchronous rotation. To obtain a more 
precise estimate of these spacings (or to see if more periodicities 
are present than suspected) will require a much longer run on this 
star. 

On the other hand, the expected spacing for consecutive m 
modes in an i = 1 multiplet is 13.60 ± 0.34 /iHz, where the uncer- 
tainty comes from the range of Cn,i for relevant models. Table |2] 
shows two spacings that are consistent with I = 1 rotational split- 
ting at the orbital frequency. The Am = 2,1 = 1 spacing involves 
the closely spaced pair fl4,15 on the high frequency side. Com- 
plicating this picture further is the observation that the separation 
between fl5 and fl2 is (to well within measurement error) pre- 
cisely equal to the orbital frequency. This suggests some geometric 
component to the overall frequency structure. 

With many periodicities present in a narrow frequency 
range, the question arises as to whether the number of pulsation 



Table 2. Frequency spacings in K1C02991403. 



IDs Spacing [/iHz] Comment 



Exact orbital frequency spacings 
fl5-fl2 26.122 (0.031) /orb - 0.002 /tHz 
fl3-fll 26.046(0.043) /orb - 0.079 aiHz 

close spacings near run resolution 
f3 - f2 0.469 (0.042) 
f5 - f4 0.548 (0.037) 
f9-f8 1.160(0.028) =2 x 0.580 

fl2-fll 1.080(0.033) =2 x 0.540 

closely spaced triplet 



fl4 


-fl3 


0.585 (0.033) 






fl5 


-fl4 


0.572 (0.029) 








I = 


1 splitting = 13.60 ± 


0.34 /tHz 




fl4 


-fll 


26.630 (0.031) 


Am = 


2 


fl5 


-fll 


27.202 (0.038) 


Am = 


2 


f8- 


f7 


13.410 (0.038) 


Am = 


1 




I = 


2 splitting = 21.91 ± 


O.lO/tHz 




flO 


-f4 


65.750 (0.036) 


Am = 


3 


fll 


-f6 


43.751 (0.030) 


Am = 


2 




I = 


3 splitting = 24.06 ± 


0.10/tHz 




fl2 


-f9 


23.959 (0.025) 


Am = 


1 


fll 


-f8 


24.039 (0.034) 


Am = 


1 



modes available to the star is sufficient to explain the number 
of periodicities seen. However, unlike the short-period sdB pul- 
sators, this is no t an issue for IpsdBV pulsators (e.g. Fig. 7 in 
iReed et alj|2004h . S tel lar models show many closely spaced fre- 
quencies in this region dFontaine et al .l2006l : |jefferv" and Sai 0I2OO7I : 
Ivan Grootel et al.ll201(J) . Still, pulsation models have found it eas- 
ier to drive high-degree modes and since the detected amplitudes 
continue down to the detection limit, and there are likely pulsations 
below our limit, it may be that high-degree modes are present. Ke- 
pler can test this critical model assumption in that with more data 
and a lower detection threshold, it may be possible to determine if 
/ > 3 are ind eed present. As Z ^ 3 have a large degree of geometric 
cancellation dCharpinet et al.ll2005l : iReed et al.ll2005h . if their am- 
plitudes are intrinsically similar to low-degree modes, then their 
observed amplitudes must be significantly reduced. 



3.3 KIC 11179657 

As for KIC 02991403, we determined the 4-a detection limit for 
KIC 11179657 using apparently variation-free frequency regions 
on either side of the range where obvious periodicities are present, 
and interpolating in between. On the low frequency side, near 
100 ^Hz, the 4-a level is 0.466 mma, and the noise level falls 
slightly to 0.431 mma near 400 /jHz. The blue lines in Fig [3] show 
the detection limit and the prewhitening process for KIC 1 1 179657. 
Vertical arrows denote periodicities removed from the temporal 
spectrum at each stage. The lowest-frequency variation has a fre- 
quency near 29.4 ^Hz and a period near 0.39 d. This frequency and 
its associated harmonic are attributed to binarity and shown in the 
left panels of Fig |3] We detected 8 other frequencies with confi- 
dence, and as in the previous section we attribute these to long- 
period (7— mode pulsations. The identified periodicities are listed in 
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Figure 2. Temporal spectra of KIC 02991403. The left panels show the peaks attributed to binary effects. The top panel 
shows the temporal spectrum of the data with the clear signature of the orbit. The middle panel has had that principal 
periodicity removed, exposing the first harmonic of the orbital period. The bottom panel shows the small residual spectral 
features after removal of the fundamental and first harmonic. The panels on the right show the peaks attributed to long- 
period pulsations in the sdB star. Again, the top panel shows the temporal spectrum of the observed light curve; the 
middle panel shows the residuals after removal of the highest eight peaks, and the bottom panel shows residuals following 
removal of all peaks in Table[T] The blue (sloping horizontal) line denotes the A-a detection threshold. 



Table[3l Prewhitening effectively removes nearly all power near the 
detection limit. 

Two peaks in the bottom panel of Fig. [3] indicate other peri- 
odicities just below the 4-(t detection limit, but likely intrinsic to 
the star. We list these as tentative peaks in Table [3] but for this pa- 
per we adhere to the 4-(j limit in identifying secure periodicities. 
There may indeed be other periodicities present, but we defer iden- 
tification of further peaks to upcoming studies with longer time- 
series data, which we anticipate obtaining in the future with Kepler. 
Again, this is likely as Kepler continues to monitor this star; longer 
runs may allow us to detect high-degree modes which suffer large 
geometric cancellation. 

In Table|4]we list spacings that correspond to multiples of the 
orbital frequency, close splittings, and spacings expected for syn- 
chronous rotation. A single splitting near the orbital frequency is 
detected. This is between f6 and f8, with a spacing of 29.66 /xHz, 
which is, within the formal frequency resolution of the data, equal 
to the orbital frequency. Among the remaining periodicities, there 
is the suggestion of a closely spaced triplet (s2, f3, f4) with an aver- 
age spacing of 0.775 /iHz, with one of the three components being 
the tentative peak s2. If this is indeed a rotationally split triplet, 
then KIC 1 1 179657 rotates with a period of 7.62 d (if / = 1). This 
is long compared to the orbital period. 

We see no spacings in the current set of frequencies that could 
be an i = 1 rotationally split multiplet at the synchronous ro- 
tation period. There are frequency differences that are at the ex- 
pected spacing for Z = 2 multiplets. Interestingly, two of the 



Table 3. Frequencies, periods, and average amplitudes for KIC 1 1 179657. 
The binary period and harinonic are included as well. Forinal least-squares 
eiTors are in parentheses. 



ID 


Frequency [/^Hz] 


Period [s] 


Amplitude [mma] 




Binary Period and first harmonic 


forb 


29.3403 (0.0024) 


34082.8 (2.8) 


9.74 (0.10) 


^forb 


58.671 (0.032) 


17044.3 (9.4) 


0.74 (0.10) 




Pulsation Frequencies 




fi 


186.517 (0.018) 


5361.43 (0.53) 


1.115 (0.088) 


f2 


195.700 (0.016) 


5109.87 (0.42) 


1.275 (0.086) 


f3 


284.656 (0.025) 


3513.01 (0.31) 


0.864 (0.088) 


f4 


285.393 (0.034) 


3503.94 (0.42) 


0.626 (0.086) 


f5 


295.690 (0.044) 


3381.92 (0.50) 


0.468 (0.086) 


f6 


308.597 (0.030) 


3240.47 (0.31) 


0.695 (0.086) 


f7 


337.167 (0.014) 


2965.89 (0.13) 


1.446 (0.086) 


f8 


338.256 (0.028) 


2956.34 (0.24) 


0.755 (0.086) 




Additional periodicities near 4-(t 


cutoff 


si 


231.634 (0.037) 


4317.15 (0.70) 


0.447 (0.086) 


s2 


283.843 (0.044) 


3523.07 (0.54) 


0.451 (0.086) 



closely spaced frequencies, f3 and f4, are (within the errors) at 
4 times the expected / — 2 spacing from fl. Again, though, like 
KIC 02991403, the high frequency component is split by a smaller 
value (the close spacing noted above). Considering these complica- 
tions and the limits of signal-to-noise and frequency resolution, no 
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Figure 3. Same as Figure[2]but for KIC 1 1 179657. Arrows indicate frequencies of those periodicities tliat were removed 
at each stage of prewhitening. The frequencies are listed in Table[3] The horizontal blue line denotes the 4-cr detection 
threshold. 



Table 4. Frequency spacings in KIC 1 1 179657. 



IDs 



spacing [/^Hz] 



Comment 



Exact orbital frequency spacings 
f8-f6 29.659(0.041) /orb + 0.319 AtHz 

close spacings near run resolution 
f3-s2 0.813 (0.050) 
f4 - f3 0.737 (0.042) 
f8-f7 1.087 (0.031) 

I = 1 splitting = 15.27 ± 0.38 ^tHz 



I : 
f6-s2 
f3 -fl 
f4-fl 

I ■■ 
f4-sl 
f8 -f2 



: 2 splitting = 24.61 ± 0.12 ^Hz 
24.754 (0.053) Am = 1 

98.14(0.024) Am = 4 

98.88 (0.038) Am = 4 

: 3 splitting = 27.05 ± 0.12 ^tHz 
53.76 (0.050) Am = 2 

54.41 (0.052) Am = 2 



firm conclusions about rotation in KIC 11179657 can be drawn at 
this point. 



4 DISCUSSION 

4.1 Synchronous rotation and rotational splitting 

The binary nature of these two stars, and in particular the short or- 
bital periods, promises to provide the extra information necessary 
for identifying rotationally split multiplets under the assumption 
that the primary star - the sdB - is forced into synchronous rotation. 
If that is the case, the discussion in Section 2 leads to the expecta- 
tion that rotationally split multiplets should be seen with splittings 
of 13.6(21.9) /iHz for I = 1(1 = 2) modes in KIC 02991403. For 
KIC 11179657, the respective expectations are 15.3 jjHz for / — 1 
modes, and 24.6 /jHz for I — 2 modes. 

Both stars show fine structure (that is, frequency separations 
that are small compared to the orbital frequency) that resembles ro- 
tational splitting for slow rotation. Identification of almost-equally 
spaced triplets in these stars suggests that these splittings are "per- 
manent" rather than a manifestation of slow amplitude or phase 
modulation. However, both stars also show larger splittings at the 
expected values for a synchronous rotation - for / = 1 in the case 
of KIC 02991403, and for Z = 2 in KIC 1 179657. In both cases, the 
spacings imply that we are seeing sectoral I = ±m and, curiously, 
the higher frequency component shows fine structure (as defined 
above). This splitting hierarchy has been seen in other g-mode pul- 
sators - in part icular, the pulsating hot white dwarf (GW Vir) star 
PG2131-F066 i lKawaleretai]|l995h . 

The relatively small number of periodicities identified (so far) 
in KIC 02991403 and KIC 11179657 makes it difficuh to identify 
multiplets with splittings compatible with the binary frequency, es- 
pecially considering that the splitting depends on the value of Cnt ■ 
We do hope that further observations will reveal more periodicities. 
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and expect that theoretical models of this star can allow us to make 
more precise estimates of Cni ■ 

4.2 Asymptotics and mode identifications 

When analysing g-modes, we often take advantage of asymptotic 
behavior of the periods to try to estimate the values of I and n. In 
particular, for p-modes, in the absence of strong mode trapping, the 
periods of modes with the same I and consecutive n are equally 
spaced in period, with a fundamental period spacing Ho that is de- 
termined by the global structure of the star. In particular, we have: 

n„i=nx^^^. (3) 

Thus, a sequence of g-modes with consecutive values of n and 
the same value of I will be spaced equally in period. Small depar- 
tures from equal spacing result from abrupt changes (with depth) 
in stars (e.g. density discontinuities that accompany composition 
gradients), and so the spacings can be a powerful probe of the inte- 
rior structure of st ars. We see this behavior ve ry clearly in pulsating 
white dwarfs (e.g. Kawaler & Bradlevlll994h . The period spacings 
depend on / in a predictable way: modes with the same value of n 
but with 1 = 1 and 1 = 2 will have periods in a ratio of approx- 
imately y/i. In those (rare) cases where sequences of I = 1 and 
I = 2 g-modes are confirmed (i.e. through rotationally split multi- 
plets), we also see that the spaci ng ratio is indeed close to 1.7 (e.g. 
PG 1159-035. ICostaetalJl2008l) . 

Up until now these properties of g-modes have not been ap- 
plied to the long-period sdB stars because of the paucity of ob- 
served modes and/or secure / determinations, and uncertainties in 
the rotational splitting (and identification of m = periods), using 
ground-based data. In Paper III we reported on the pulsations of 
other long-period sdB pulsators that are Kepler survey targets. In 
at least one case, KIC 10670103, we do see sequences of equally 
spaced periods, and in fact can identify an Z = 1 sequence with 
a period spacing of 251 s, and wl = 2 sequence with a spacing 
of 146 s. The ratio is 1.72, and confirms the expectations from the 
asymptotic behavior of g-modes. 

For the two stars in this paper, are there periods that are 
equally spaced? Tables[5]and|6]address this question by listing the 
periods in each star, along with the period differences between con- 
secutive (identified) periodicities. Though the number of periodic- 
ities is smaller than for KIC 10670103, the periods seem to show 
similar spacings to that star The mean spacing of about 136 s in 
f4 through fl6 in KIC 02991403 is similar to the / = 2 spacings 
in KIC 10670103. This spacing is inconsistent with the preliminary 
results based on rotational splitting, which favors an Z = 1 assign- 
ment for fU and fl4 (and fl5). In KIC 11179657, with fewer pe- 
riodicities, the period structure also shows spacings that are com- 
patible with those seen in KIC 02991403 and KIC 10670103. In 
summary, period spacings are insufficiently well defined to assign 
secure values of I for these stars. This may be because rotational 
splitting is present and we are not seeing m = modes. 

Even though we cannot say, with confidence, that we see an 
equal period spacing in either of these stars, we can still use the 
asymptotics of p-modes to say something about the value of If 
the driving mechanism acts on a limited number of overtones (n), 
we would expect to see groups of periods in the temporal spectrum 
with the centroids of those groups with a period ratio of 1.73. In 
KIC 1117657, the strongest mode (f7) has a period of 2966 s. If 
that is an / = 2 mode, the corresponding I = 1 mode should have 
a period of roughly 5140 s. This is fairly close to f2 in that star. 



iBDic J. r cnoa 


spacings in 




ID 


Period [s] 


Difference [s] 


fl 


6365.35 








~ 1234.8 (= 5 X 246.95) 


(f2+f3)/2 


5130.59 




(f4+f5)/2 


4332.13 








~ 550.51 (= 4 X 137.63) 


f6 


3781.62 








269.50 (= 2 X 134.75) (f7 in between) 


(f8+f9)/2 


3512.12 








137.12 


llU 


33 /4.93 








135.82 


(fll+fl2)/2 


3239.11 








252.41 (=2 X 126.21) 


fl4 


2986.70 








276.77 (= 2 X 138.39) 


fl6 


2709.93 




iable o. reiioa 


spacings in 


JvlU ill /yoj /. 


ID 


Period [s] 


Difference [s] 


fl 


5361.43 








251.56 


f2 


5109.87 




si 


4317.15 








804.14 (=6 X 134.02) 


s2,f3,f4 


3513.01 








~ 131.09 


f5 


3381.92 








141.45 


f6 


3240.47 








~ 279.4 (=2 X 139.7) 


f7,f8 2956.34, 2965 


.9 



which is the second strongest pulsation: the ratio of f2 to f7 is 1 .723. 
This suggests that we are seeing I = 2 and I = 1 pulsations in 
this star. We note that Kepler data on a single long-period sdBV 
star, KIC 10670103 (Paper III), shows the same pattern of apparent 
/ = 1 and / = 2 period ratios. 

Similarly, among the main periodicities in KIC 02991403, the 
highest-amplitude mode has a period of 2987 s. Assigning that 
mode to / = 2 suggests that the I = 1 counterpart should be at 
a period of 5173 s, which is very close to the doublet (f2,f3) at 
5130 s. The ratio of fl4 to f2 is 1.720. This too is remarkably close 
to the asymptotic value of 1.732. This simplistic preliminary anal- 
ysis suggests that we are indeed seeing I = 2 and I = 1 modes. 
This, in turn, means that the largest amplitude modes have 1 = 2 
in some cases - which is surprising since the higher-degree modes 
suffer more geometric cancelation across the surface. 

Thus there are several lines of evidence, all circumstantial at 
this point, that we are seeing / = 1 and / = 2 g-modes in these 
two stars. We expect to have longer runs (3 months or more) on 
these stars as the Kepler Mission continues, and expect that the in- 
crease in frequency precision and signal-to-noise should help us 
determine, with certainty, the rotational state of these stars. Further 
exploration of the binary properties (i.e. a radial velocity study ver- 
ifying the orbital and rotational signatures) should provide clues as 
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to the properties of the (presumably dM) companions, and the na- 
ture of the binary light curves. With this information in hand, these 
stars should play an important role in improving our understand- 
ing of the binary nature of sdB stars in general, and the role that 
binarity plays in their formation and subsequent evolution. 
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